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A 6 term, which couples to topological charge, is added to the two-dimensional lattice CP 3 model and U(l) 
gauge theory. Monte Carlo simulations are performed and compared to strong-coupling character expansions. In 
certain instances, a flattening behavior occurs in the free-energy at sufficiently large 9, but the effect is an artifact 
of the simulation methods. 



Following the discovery of instanton solu- 
tions in four dimensional Yang-Mills theories 
pi, the importance of adding a 9 term Sg — 
jp9 J d 4 xF^F^(x)/(32ir 2 ) to the action was re- 
alized @d. 

Since So breaks parity, time-reversal 
invariance and CP symmetry when S / or 
9 ^ 7T, the strong interactions explicitly violate 
these symmetries for < < it. The physi- 
cally effective 9 angle is bounded experimentally 



by 9 eff < 10" 



The question of how 9 e ff 



can naturally be so small constitutes the strong 
CP problem in QCD. 

Due to the complexity of the problem a pre- 
liminary study of simpler systems on the lattice 
is useful. A class of such systems are the two di- 
mensional CP^ 1 models HQ], which have many 
features in common with four dimensional Yang- 
Mills theory. 

Let us start with a general analysis of sim- 
ulating systems with 9 terms. For the lattice 
U(l) gauge theory and the CP N ~ X model, the 
local topological density v v is defined via v v = 
log (Up) /(27r), where U p is the product of the 
U(l) link phases around the plaquette p and 
where — ir < log (U„) < ir. The total topologi- 
cal charge Q is given by Q — J2 P v v The theta 
term Sg term is i9Q, that is, [p| 



Sg term = — ^2 l0g (^p) 



(1) 



Eq. (|l|) is the lattice analog of the continuum 9- 
term action J cPxFqi. 

Let f(9) be the difference between the free en- 



ergy J- (9) of a system with a 9 term and the free 
energy of a system with 9 = 0: 



f(9) = T{9) - ^(0) 



(2) 



Typically, f(9) is an increasing function of 9 for 
< 9 < it. For a fixed volume V, let P(Q) be the 
probability of having a configuration with topo- 
logical charge Q in the system. The free energy 
difference f(9) is then constructed from P(Q) us- 



ing 



exp (-Vf(9)) = P (^) ex P 



(3) 



Normally P(-Q) = P{Q), so that f{-9) = f(6). 

In a Monte Carlo simulation, an approximation 
/mc(#) to f(9) is obtained by using a measured 
Pmc{Q) m lieu of P(Q). Hence 

- Vf M c(6)) = log [exp (-Vf(9)) + 6Z(9)} (4) 

where SZ(9) = J2q $P{Q) exp (z6»Q).Q Since f(9) 
is an increasing function of 9, an accurate mea- 
surement of f(9) for < 9 < 9b is obtained if 



\SZ(9)\ «exp(-U/(M) 



(5) 



In particular, since /(0) = and \5Z(9)\ <C 1, 
there is always a region near 9 = for which 
f(9) can be measured in a Monte Carlo simula- 
tion. However, away from 9 = 0, Eq. (|^) implies 
that for sufficiently large V, a limiting value of 
9b exists beyond which it is impossible to reliable 



lr The deviation between Monte Carlo measurements and 
exact results is denoted by S. 
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Figure 1. U(l) Free Energy Versus 6 at (3 = 1.0 
for the Naive and Binning Methods. 



compute f(0). The value of 0s depends on the 
statistical accuracy of the simulation. As V gets 
larger, 9 b decreases unless enormous numbers of 
measurements are undertaken to reduce statisti- 
cal errors. For large V, obtaining enough mea- 
surements becomes, in any practical sense, im- 
possible. Clearly, it is more difficult to measure 
f(9) throughout the entire fundamental region of 
9, as V gets larger. 

It turns out jllj that in most Monte Carlo sim- 
ulations, there is a tendency for 

|£P(0)| > |£P(1)| > |5P(2)| > ... . (6) 

Now if \SP(0)\ is much larger than the other 
\SP(Q)\ then, from Eq. (g), one deduces an es- 
timate for 9b 

/(0 B )«A|log|M>(O)|| . (7) 
Since Monte Carlo results are reliable for 9 < 9b, 
fMc(0)~f(9) ioY9<9 B • (8) 
If, in addition, <5P(0) > 0, then one finds 

f M c(0) « ~ log SP(0) for 9 > 9 B , (9) 

so that a constant "flat" behavior in /mc($) will 
be observed, a pure artifact of the simulation. If, 
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Figure 2. CP 3 Free Energy Versus 9 at j3 = 0.2 
on 4 x 4 and 6x6 lattices. 



on the other hand, SP(0) < 0, then the measured 
/mc(S) will blow up for 9 > 9b- In our simula- 
tions we have observed both types of behaviours. 

From the above discussions we see that as long 
as finite-size effects are under control, that is 
£ < V^ 1 /^,^] small- volume results for the mea- 
surement of f{9) are more reliable than large- 
volume results. If a flattening behavior of the free 
energy f{9) for large 9 is observed, one should be 
cautious that the result is spurious. In particu- 
lar, one should try to see whether |<5P(0)| is bigger 
than the other \SP(Q)\. Therefore the guideline 
emerges that if a large- volume simulation shows 
a flattening effect for f(9) for 9 sufficiently large, 
but a smaller- volume simulation does not, one 
should trust the smaller- volume result. 

We note that in the work of || a flat behaviour 
of the free energy was observed and attributed to 
a phase transition. The results of this work fl0]| 
suggest that the flattening is a simulation effect. 

The 2-D lattice U(l) gauge theory serves as 
an ideal testing ground, as computer simulations 
can be compared to exact analytic results JljJ ]. 
Figure 1 plots the free energy versus 9 for /? = 
1.0 on a periodic 16 x 16 lattice for two different 
runs. The solid line is the exact analytic result. 

2 Here, £ is the correlation length and d is the number of 
dimensions of the system. 
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Figure 3. CP 3 Free Energy Versus at (3 = 0.6 
on 6 x 6 and 8x8 Lattices. 



Both runs have comparable statistics and agree 
with the analytic results for 9 less than 2.1, the 
value of the "barrier f?" 9b- Using the known 
error <5P(0) in Eq.(0) to estimate 9 b, one finds 
9b ~ 2.05, confirming the above data analysis. 
The run exhibiting the anomalous flat behaviour 
in the free energy for 9 > 2.1 in fact has a positive 
SP(0), as predicted by Eq. (§). 

For the simulations of the lattice CP 3 model^ 
we have employed the "auxiliary U(l) field" for- 
mulation |l3]] . Figures 2,3 and 4 show the free 
energy for (3 = 0.2,0.6 and 0.7 on 4 2 ,6 2 and 8 2 
lattices. The solid line represents the tenth-order 
strong-coupling character expansion of ref. p2| . 
Figures 2 and 3 show that simulations on smaller 
lattices are more reliable, as the simulations on 
the larger lattices exhibit anomalous flattening. 
Again the estimated 9b for these simulations was 
in good agreement with the observed one. In the 
intermediate coupling regime of (3 — 0.7 in fig- 
ure 4 the Monte Carlo data is most likely to be 
trusted over the strong-coupling expansion. Cu- 
riously for higher values of (3 the MC simulations 
were nicely fitted by a cosine JTo[ ] , which also 
arises from a topological gas picture 




Figure 4. CP 3 Free Energy Versus 6 at /? = 0.7 
on a 8 x 8 Lattice. 
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3 For simulations without a 8 term see refs. in fli 



